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ABSTRACT 

We present results from a 42 ks Chandra / ACIS-S observation of the transitional FRI/FRII 
radio galaxy 3C 288 at z = 0.246. We detect ~3 keV gas extending to a radius of ^0.5 Mpc with 
a 0.5-2.0 keV luminosity of 6.6 x 10 43 ergs s , implying that 3C 288 lies at the center of a poor 
cluster. We find multiple surface brightness discontinuities in the gas indicative of either a shock 
driven by the inflation of the radio lobes or a recent merger event. The temperature across the 
discontinuities is roughly constant with no signature of a cool core, thus disfavoring either the 
merger cold-front or sloshing scenarios. We argue therefore that the discontinuities are shocks 
due to the supersonic inflation of the radio lobes. If they are shocks, the energy of the outburst 
is ^10 60 ergs, or roughly 30% of the thermal energy of the gas within the radius of the shock, 
assuming that the shocks are part of a front produced by a single outburst. The cooling time of 
the gas is ~10 8 yrs, so that the energy deposited by the nuclear outburst could have reheated 
and efficiently disrupted a cool core. 
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1. Introduction 

Chandra has revolutionized our knowledge of 
the energetic processes involved in the jets, the 
lobes and the nuclei of radio galaxies. The X- 
ray images of the hot atmospheres in galax- 
ies show a wealth of structures associated with 
central radio sources, including cavities, metal 
enriched plumes, filaments, and shock fronts 



(IPeterson fc Fabian 



2007 



200 6t iMcNamara fc Nulsen 
McNamara et al. I I2009F ) 



High-resolution 
spectroscopic observations from Chandra and 
XMM-Newton have conclusively ruled out simple, 
steady cooling flow models a nd this has been one 



of the significant discov ery (|Peterson fc Fabian 
20061 : iDavid et aL~ll2006l) . Since the cooling time 
of gas in many cluster cores is much less than the 



Hubble time, energy must be occasionally or con- 
tinuously supplied to cluster cores to prevent the 
formation of cooling flows. A primary candidate 
for the suppression of cluster cooling flows is feed- 
back between the central supermassive black hole 
(SMBH) of active galaxy and the cluster gas. 

Studies of the X-ray gas environments of FRI 
and FRII sources with Chandra have led to im- 
portant new constraints on how jets propagate 
through their ambient media, and how radio lobes 
interact with and transfer ener gy to their large- 
scale gas environments (e.g ., Hardcastle et al. 

Kraft et al. I l2006h 



2002; Sambruna et al. 200- 



However, 'high-excitation' FRII sources tend to 



shifts (lEllinESon et al. Ill991a 


b c : Harvanek et al. 


2001. 20021 Isobe et al. 120051: 


Kraft et al. 112007b). 
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unless the mechanical p ower of the jet is unusu- 
ally high (e.g., Cygnus A lSmith et aUhOO'A . This 
strongly suggests that the jet power and the in- 
teraction between jets and the hot gas in their 
vicinity play pivotal roles in governing the overall 
morphology (FRI vs. FRII) of the radio source. 

The radio source 3C 288 is identified with an 
elliptical galaxy with m v = 18.3 (M v — —22.08, 
Goodson et al. I 1979() . It is the archetypal exam- 



in Section 4. Section 5 contains a brief sum- 
mary and conclusions. We assume WMAP cos- 



pie of a trans itional FRI/FRII ('jett ed double') 
radio galaxy (jFanaroff fc Rilev I Il974j ). In opti- 
cal imaging, there are fainter galaxies in the field 
which are presumed to be me mbers of a cluste r 
that is dominated by 3C288 (jWvndham I Il966l ). 
Dominant cluster members can produce asym- 
metric "wide-angle- tail" (WAT) radio structures 
(|Burns et al. | ["l986). despite their presumed low 
peculiar velocities, but 3C288 is not a clear ex- 
ample of WAT morphology Its monochromatic 
power at 1.5 GHz is 2.5 x 10 33 ergs s -1 Hz -1 
and its integrated spectral in dex an'^ between 



0.75 GHz and 5.0 GHz is -0.97 (|Laing fc Peacock 



198(1 s v oc v a ). Its radio luminosity would place 
it firmly in the high-power (i.e., FRII) regime, 
and yet its observational characteristics, in partic- 
ular its edge-darkened radio morphology, are strik- 
ingly different from those of the canonical 'clas- 
sical double' radio galaxies. Although, the edge 
darkening of the radio structure, and its spec- 
tral distribution, are reminiscent of a (distorted) 
radio "trail", no other properties of 3C288 sug- 
gest that it belongs to this morphological class. 
The structure of 3C 288 is more asymmetric than 
those of nor mal double sources of its size and 
radio power ( Bridle et al. 19891 ). VLA observa- 



tions reveal a jet and a counterjet near the ra- 
dio core, and faint "wings^ of emission connected 



to the elongated lobes (jBridle et al. I Il989h . The 



small radio size, unusual morphology, and polar- 
ization asymmetries of 3C 288 raise the possibility 
that it is interacting s trongly with the ambient 



gas ((Bridle et al. Ill989h . Thus, 3C288 is an ideal 



candidate for study with Chandra to better un- 
derstand the role that AGN heating may play in 
the formation of radio structure in the hot phase 
of the intergalactic medium (IGM). 

This paper is organized as follows: Section 2 
contains a summary of the observational de- 
tails. Results of the data analysis are presented 
in Section 3 and we discuss their implications 



mology throughout this paper ([Spergel et al. 
20071 ). The observed redshift (z = 0.246) of 



the host galaxy of 3C 288 corresponds to a lu- 
minosity distance of 1192.9 Mpc, and 1" = 3.725 
kpc. All coordinates are J2000. The elemental 
abundance s are relative to the Solar value tab- 
ulated by lAnders fc Grevassel ( 19891 ). Absorp- 
tion by gas in our galaxy (Nh = 9 x 10 19 cm~ 2 , 
Dickey fc Lockman Il990h is included in all our 
spectral fits. 

2. Observations 

The radio galaxy 3C 288 was observed on 2008- 
04-13 (OBSID 9275; P.I. D.A. Evans) with Chan- 
dra/ ACIS-S in VFAINT mode for -42.0 ks. We 
made light curves for each CCD in the 0.5 keV to 
10.0 keV band in order to search for background 
flares and intervals where the background rate 
was high; none were present, leaving 39647.8 s 
of good data. We performed the usual filtering 
by grade, excluded bad/hot pixels and columns, 
removed cosmic ray "afterglows" , and applied 
the VF mode filtering using tools built in CIAO 
(http ://cxc.harvard.edu/ciao] ). Images were gen- 
erated after subtracting background and correct- 
ing for exposure (which included all the effects 
mentioned above). 

We use archival VLA observ ations of 3C 288 at 
4.885 GHz (jBridle et al. Ill989h - The map is taken 



from the online 3CRR AtlaH which provides well- 
calibrated, well-sampled images. 

3. Data Analysis 

A Gaussian-smoothed (radius fwhm = 2") 
Chandra/ ACIS-S image of 3C288 in the 0.23- 
5.00 keV band with radio contours overlaid is 
shown in Figure [TJ All point sources other than 
the 3C 288 nucleus have been removed. We detect 
diffuse thermal emission from the cluster ICM 
with a temperature of —3 keV extending to —390 
kpc (~104"). Diffuse emission from the source fills 
only a fraction (—3%) of the area in the S3-chip, 



2 ATLAS catalog: Radio images and other data for the near- 
est 85 DRAGNs (radio gala xies and related obj ects) in the 
so-called "3CRR" sample of lLaing et al. (1983). Available 
at http://www.jb.man.ac.uk/atlas/index.html . 
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so local background was extracted for all spectral 
analyses, from the source-free region on the S3 
chip. The radio source is small, ~ 67.1 kpc (= 
18" across), and the X-ray bright gas core lies ~ 
11.2 kpc (~3") north of the radio core; this off- 
set is perhaps associated with the non-hydrostatic 
motion of the gas. The X-ray isophotes within 
112 kpc (0.5') of the nucleus are circular. For 
simplicity we assume spherical symmetry in our 
analysis below. On larger scales, however, the 
isophotes show an extension to the south-east. 

Interestingly, we detect significant X-ray emis- 
sion below 0.50 keV. Such emission is not com- 
monly observed in clusters of galaxies with Chan- 
dra, and was initially suggestive of inverse Comp- 
ton scattering of CMB photons from a large 
population of cosmic ray electrons in the clus- 
ter core , similar to that claimed for the Coma 
cluster ( Sarazin fc Lieu 1998t Finoguenov et al. 



120031: lErlund et al. 1120071 ). This emission lies pre- 
dominantly at the cluster core. Additionally, it 
appears that there is a cavity in the gas asso- 
ciated with the southern radio lobe and that the 
0.23-0.50 keV emission is roughly aligned with the 
jet axis. However there is no obvious direct cor- 
respondence between this soft X-ra y emission and 



radio features as seen in Hydra A (jNulsen et al. 



2005|). An image of the cluster in the 0.23-0.50 
keV band with radio contours overlaid is shown in 
Figure [2l We created surface brightness profiles 
of the emission in the 0.3-0.5 keV (soft) and 0.5- 
1.5 keV (hard) bands, subtracted an appropriate 
background, and fitted a line to the ratio of soft 
and hard X-ray emission as a function of distance 
from the center. The slope of this line is con- 
sistent with zero within the uncertainties and we 
confirm that, within the statistical uncertainties, 
there is no difference in the spatial distribution 
of this soft X-ray emission relative to the hard 
X-ray emission. Spectral analysis of the Chandra 
data, combined with spectral analysis of archival 
ROSAT PSPC data, confirm that this emission 
is simply the bremsstrahlung continuum emission 
from a ~3 keV plasma. It is only visible in the 
3C 288 cluster because of the combination of gas 
temperature and unusually low Galactic column. 

3.1. Compact Components 

Contours from the 4.9 GHz radio map of 3C 288 
are shown in Figure [T] overlaid on the background- 
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Fig. 1. — Gaussian-smoothed (radius fwhm = 
2"), exposure corrected, background subtracted 
C7iandra/ACIS-S image of 3C288 in the 0.5-2.5 
keV band. All point sources, other than the ac- 
tive nucleus of 3C288, have been removed. We 
see diffuse thermal emission with a temperature of 
~3 keV extending to ^510 kpc (~1.73'). Contours 
from the 4.9 GHz radio map of 3C 288 are overlaid 
and 10 contour levels are placed linearly between 
0.4 and 12.0 mJy beam -1 . The radio source is 
small, 294 kpc (1.3' across) and the radio core is 
conincident with the optical host galaxy. A peak 
of X-ray emission lies ~11.2 kpc (~3.0") east of 
the optical host galaxy. 





















Fig. 2.— Chandra/ ACIS-S image of 3C 288 in the 
0.23-0.50 keV soft X-ray band with 5 GHz radio 
contours (0.6" resolution) overlaid. This image 
shows the soft X-ray excess, which is not com- 
monly observed in clusters of galaxies and is uni- 
formly detected across the whole 42 ks observa- 
tion. 
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subtracted, exposure-corrected, Gaussian-smoothed 
Chandra ACIS-S image. The core is the brightest , 



most compact radio feature ( Bridle et al. 1989h 



and c oincides with the opti cal position of the nu- 



cleus ( Goodson et al. Ill979l ). It has a flux density, 
5 4 .9GHz, of 12.68 ±0.07 mJy (1 mJy = 10~ 26 
ergs s _1 cm~ 2 Hz" 1 ). The northern and south- 
ern hotspots are clearly detected in 4.9 GHz and 
1.4 GHz radio maps and have flux densities of 6.88 
±0.09 mJy and 5.84 ±0.11 mJy, respectively at 
4.9 GHz. The radio core has a spectral index of 
a iA §Hz = -°- 76 ±0 - 02 i s u °c v a , where S v is 
the flux density at frequency V) and the error bars 
are la confidence). The northern and southern 
hotspots have spectral indices of af;° = —0.76 
±0.04 and -0.87 ±0.06, respectively. We do not 
detect an X-ray point source (see Figure [J) at the 
location of the radio core. Only eight net counts 
above background were extracted in a 0.5" radius 
circle at the location of the core with no point 
source evident above the level of the diffuse emis- 
sion from the gas. The 3cr upper limit to the 1 keV 
(rest frame) flux density of the core is ~0.31 nJy 
assuming a photon index of 2 and Galactic absorp- 
tion. Assuming that this upper limit represents 
nonthermal X-rays from the AGN, the upper limit 
to the radio to X-ray power-law index is OL T ~^ y ~ 
—0.99 ±0.02. The corresponding core radio lumi- 
nosity at 178 MHz assuming a flat spectrum and 
the unabsorbed core X-ray luminosity at 1 keV 
are consistent with the expectation from the cor- 
relation between these two quantiti es shown in 
Hardcastle. Evans, fc Croston I (l2009h . 

No X-ray emission is detected from the jet or 
the compact hot spots of the two lobes. With 
no point source evident and zero counts above 
the diffuse emission from the gas, we use the 3c 
counts of seven, by measuring the off-source back- 
ground level and then use Poisson statistics to 
find the number of counts corresponding to a 3- 
sigma Gaussian probability, for the northern and 
the southern hotspots. The 3er upper limit to the 
1 keV (rest frame) flux density of both, the north- 
ern and southern hotspots is ^0.29 nJy assuming 
a photon index of 2 and Galactic absorption. We 
deduce upper limits to the X-ray-radio power-law 
indices for the northern and southern hotspots of 

«fadio ay > -°- 96 ±0 - 02 and > -°- 95 ±0 - 02 > respec- 
tively (again the error bars are at ler confidence). 
If we use these spectral indices, the expected flux 



densities in the optical and infrared bands, S ° 

5000 A 

and 5*7^ are 0.09 and 1.15 fiJy for the northern 
hotspot, and 0.09 and 1.07 /LtJy for the southern 
hotspot, respectively. These are undetectable with 
the current generation of optical and infrared ob- 
servatories. There is no non-thermal emission de- 
tected from the nucleus, lobes, or jets from this 
radio galaxy, and the upper limits are consistent 
with detections of such emission in much closer 
radio galaxies. 

3.2. Large-Scale Diffuse X-ray Emission 

The extended, diffuse X-ray emission seen in 
Figure [T] is attributed to emission from the hot 
gas of a cluster atmosphere. We derived a global 
temperature and metallicity for 3C 288 within a 1' 
radius circular region covering the majority of the 
cluster emission. The spectrum was extracted us- 
ing the CIAC0 specextract tool, binned to 10 cts 
per bin and fitted in the 0.5-5.0 keV range using an 
abso rbed APEC m odel within the XSPEC pack- 
age ( Arnaudlll996l) . The neutral hydrogen col- 
umn density was fixed at the Galactic foreground 
value of A H = 9.0 x 10 19 cm~ 2 . The best fitting 
values for temperature and abundance are fceT 
= keV and Z = 0.73±g;i|, respectively, 

where the errors are 90% confidence limits. If we 
allow the neutral hydrogen column density to be a 
free parameter, the best fitting values for temper- 
ature and abundance do not change appreciably, 
and the changes are smaller that la uncertainties. 
The unabsorbed X-ray luminosity from the best 
fitting model in the energy range 0.5-5.0 keV is 
1.11 ±0.03 x 10 44 ergs s^ 1 within the r = 223.5 
kpc (1') circular aperture. Additionally, the unab- 
sorbed X-ray luminosity in the energy range 0.5- 
2.4 keV is 7.60 ±0.18 x 10 43 ergs s" 1 and the re- 
sult is c onsistent with the expectations for a 3 keV 
cluster (jMarkevitch 1119981 ) . The azimuthally aver- 
aged radial surface brightness profile of the X-ray 
emission from the gas is shown in Figure [3] The 
best-fitting isothermal /3-model profile has been 
overlaid. We find f3 = 0.52 ±0.02 and a core ra- 
dius ro = 11.94" ±0.93" from fitting the surface 
brightness profile between 3" and 200" from the 
nucleus. 

We find two breaks in the surface brightness 



3 All spectral extraction and spectral analsis were performed 
using CIAO v4.1, the CALDB v4, and XSPEC v!2.5.1. 
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Fig. 3. — Azimuthally averaged radial surface 
brightness profile of the X-ray emission from the 
gas in the energy band from 0.5 keV to 2.5 keV. 
The best-fit isothermal /3-model profile has been 
overlaid. The two vertical lines show the location 
of prominent, sharp edges in X-ray brightness at 
distances from the nucleus of 53.3 kpc (14.2") and 
91.1 kpc (24.4") along the eastern and the south- 
western directions, respectively. 



distribution, one 14.2" east of the nucleus and 
another 24.4" to the southwest, as shown by 
the black arrows in Figure 2J The presence of 
these breaks in the surface brightness distribu- 
tion implies a sharp change in the gas density 
or temperature of the gas across the disconti- 
nuities. Chandra has observed a large num- 
ber of simil ar features in other clus ters, such as 
Abelll795 (iMarkevitch et al. I l200ll). Abell 3667 




2( 



Vikhli nin et al. I l2001aD, M87 (IForman et al 



on the L x -T relation. Thus, we consider this 
po ssibility to be unlikely. Followi ng the analysis 
of Markevitch fc V ikhlinin ] ioo3), below we de- 
termine the temperature and pressure across the 
discontinuity to evaluate which of the other two 
scenarios, sloshing or supersonic inflation of radio 
lobes, is more plausible. 




Fig. 4. — Gaussian-smoothed (radius fwhm = 
2") Chandra/ ACIS-S image in the 0.23-2.50 keV 
band. The image has 2" bins. The black sets of 
arrows denote the positions of two surface bright- 
ness discontinuities. The image also show re- 
gions where surface brightness profiles were fitted, 
and where spectra were extracted for temperature 
measurements. A large 1' region centered on the 
radio galaxy was used for global temperature and 
abundance measurements. The surface brightness 
profiles shown in Figures 5 and 6 are taken from 
the two sectors. 



2005Ll2007h. HvdraA (iNulsen et al. 112005b. MS 0735.6+ 7421 

al 



( McNamara et al. |2005j) , and Abell 1201 (jOwers et 
2009), and they arc generally attributed to three 
phenomena: merger cold fronts, sloshing of clus- 
ter cores due to non-hydrostatic motions of the 
gas, and shocks due to nuclear outbursts. Most 
merger cold fronts are offset from the center (e.g. 
Abell 3667) and are the result of the infall of a 
massive subclump into the cluster. Unless we 
are viewing such a merger head-on, the morphol- 
ogy of these features in the 3C 288 cluster gas is 
very different than what we observe in Abell 3667 
(jVikhlinin et al. I l2001af ) . Additionally, it would 
be surprising, if we are witnessing a major merger 
from such a viewing angle, that the cluster lies 



We fitted absorbed APEC models to three an- 
nuli in two sectors (shown in Figure |4]) centered 
on 3C288. The vertex of the annuli was fixed at 
the nucleus, but the binning of the annuli was ad- 
justed, so that the radius of curvature of the sec- 
ond and third annular bins in each fit matches that 
of the associated discontinuity (i.e. we created the 
bins so that two annuli were interior and one an- 
nulus was exterior to the discontinuity). The goal 
of this spectral fitting was to determine whether 
the gas temperature interior to the discontinuities 
was hotter or cooler than the exterior gas temper- 
ature. Only the temperatures and normalizations 
were free parameters in these fits. The elemental 
abundance was frozen at the best-fit value deter- 
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mined in the global fits (Z = 0.73). Plots of the 
temperature profiles for the two different sectors, 
one between position angle (P.A.) = 30° and P.A. 
= 130° and another between P.A. = 180° and P.A. 
= 260°, as a function of radius from the phase- 
center (position of the host galaxy) are shown in 
Figures [5] and [6j panel (c) . We find no signifi- 
cant jump in the projected temperature across ei- 
ther discontinuity within 90% error uncertainties 
of AT/T -37%. 

We deprojected the surface brightness profiles 
(using the best-fit gas temperature derived above) 
to determine the density and pressure as a function 
of distance from the nucleus. Figure 0] shows the 
ACIS-S image of the central region, with promi- 
nent, sharp edges in X-ray brightness at distances 
from the nucleus of 53.3 kpc (14.2") and 91.1 kpc 
(24.4") along the eastern and the south-western 
directions, respectively. We model each of these 
surface brightness discontinuities with a broken 
power-law density model. We fitted the surface 
brightness across the discontinuity in the sectors 
shown in Figure|4]in the energy range 0.5-5.0 keV. 
The deprojected density and pressure profiles be- 
tween 30° and 130° (eastern direction), and be- 
tween 180° and 260° (south-western direction) as 
a function of distance from the cluster center are 
shown in Figures [5] and El panels (b) and (d) , re- 
spectively. Both brightness profiles have a charac- 
teristic shape corresponding to a projection of an 
abrupt, spherical (within a certain sector) jump 
in the gas density. Best-fit radial density mod- 
els of such a shape are shown in panel (b), and 
their projections are overlaid on the data in panel 
(a) of Figures [5] and El From the amplitude of 
the best fitting surface brightness model, we de- 
rived a density jump of 1.48^o 24 for the eastern 
discontinuity and 1.75^Q'gg for the south-western 
discontinuity. The confidence ranges for the den- 
sity jumps were computed from the extremes of 
the 90% confidence ranges for the best fitting sur- 
face brightness model. As confirmation, the den- 
sity profile that is derived from the isothermal, 
azimuthally symmetric surface brightness profile 
is broadly consistent with the deprojections of the 
broken power-law models. 

To minimize projection effects, we would like 
to measure the deprojected temperature profile 
across the jumps; however, the limited number of 
photons simply does not permit this, although the 



lack of temperature variation suggests that pro- 
jection effects are probably not large. Instead, we 
used the extracted spectra from two large regions 
for each front for the subsequent analysis; one for 
the bright side of the front (inside or post-shock) 
and one for the faint side of the front (outside 
or pre-shock). For the east jump, we measured 
a temperature of 2. 88^03 keV inside the jump 
and 3.461q 53 keV outside the jump, and for the 
south-western jump, we measured a temperature 
of 3.99;^3 keV inside the jump and 3.6ltJ;^ keV 
outside the jump. The temperature ratios for the 
two discontinuities are 0.83:^23 and L11 ±a4i for 
the east jump and the south-western jump, respec- 
tively. These uncertainties are at 90% confidence. 

4. Interpretation 

4.1. Are the Discontinuities Shock Fronts 
or 'Sloshing' ? 

There are at least two possible explanations for 
the surface brightness discontinuities: they could 
be shocks in the gas due to the supersonic inflation 
of the radio lobes, or they could be contact discon- 
tinuities created by non-hydrostatic motions of the 
gas core due to a recent merger (i.e. 'sloshing'). 
Unfortunately the quality of the data is not suffi- 
cient for us to make a definitive statement about 
whether these discontinuities are shocks or slosh- 
ing. The temperature of the cluster is sufficiently 
hot (—3 keV) that the cluster temperature is not 
well constrained with less than several thousand 
counts per spectrum. We favor the shock model, 
as described below, but arguments can be made for 
either scenario. We describe some of the implica- 
tions of both models. Note that our conclusions 
about the overall energetics of the radio galaxy are 
not significantly different under the two scenarios. 
More explicitly, the estimate energy of the out- 
burst is only tens of percent larger for the shock 
scenario than the subsonic inflation scenario. 

4.2. Supersonic Inflation of the Radio 
Lobes 

Several lines of evidence suggest that the dis- 
continuities are shocks, including the relative sym- 
metry of the discontinuities around the nucleus 
which are not typically seen in sloshing cores, the 
lack of cool gas in the cluster core, the lack of any 
evidence of a recent merger in the HST image of 



G 




Fig. 5. — Shock in 3C288 showing the eastern edge (sector between P.A. = 30° and P.A. = 130°). Panel (a) 
shows X-ray surface brightness profile. The filled circles are the model values that correspond to the best-fit 
gas density model shown in panel (b). Panel (c) shows the temperature profile. Panel (d) shows pressure 
profile obtained from the temperatures on either side of the edge and density profile. Error bars are 90%; 
vertical dashed line in panel (a) show the position of the density jump. 



the host galaxy (|de Koff et al. Ill996l ). the lack of 
any evidence of merging on larger scales, and the 
transformational morphology of the radio source 
from FRII to FRI suggesting that it is strongly 
interacting with its environment. If these features 
are shocks, we can estimate the expansion velocity 
of the radio lobes, total power of the outburst, etc. 
to constrain the impact the effect of this outburst 
would have on the surrounding gas. These conclu- 
sions could be definitively confirmed (or refuted) 
with a deeper Chandra or XMM-Newton observa- 
tion. 

We use the example of M87 to provide a nearby 



analogy to interpret the features seen in the 
3C288 cluster gas. M87 is the dominant cen- 
tral galaxy in the Virgo cluster showing remark- 
able morp hological details in both the X-ray and 



20001 : iHines et al. 1 1 19891 



the radio (lOwen et al. 

Forman et al. I l2005h . Radio observations show 
evidence of two distinct nuclear outbursts, and 
X-ray observations show at least two sets of sur- 
face bright ness discontinuities in the gas indicative 
of shocks ( Forman et al. 1 12007). 3C288 shows a 



signature of a shock, two regions of enhanced pres- 
sure, at a radius of 53.3 kpc on the east and 91.1 
kpc on the south-west, similar to M 87. 
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Fig. 6. — X-ray properties in the south-western edge (sector between P.A. = 180° and P.A. = 260°). Panel 
(a) shows X-ray surface brightness profile. The filled circles are the model values that correspond to the 
best- fit gas density model shown in panel (b). Panel (c) shows the temperature profile. Panel (d) shows 
pressure profile obtained from the temperatures on either side of the edge and density profile. Error bars 
are 90%; vertical dashed line in panel (a) show the position of the density jump. 



If these features are indeed shocks, the post- 
shock gas temperature and density will be higher 
than the pre-shock values in a narrow region be- 
hind the shock. With the existing data, we can- 
not determine the gas temperature profiles across 
the discontinuities with sufficient accuracy to con- 
firm that they are shocks. Panel (c) in Figures [5] 
and [6] shows the gas temperature profiles across 
the edges. For a shock discontinuity, the Rankine- 
Hugoniot jump conditions directly relate the gas 
density jump, r ~ p s hock//Oo, and the tempera- 
ture jump, t ~ T s hock/2bj where subscripts '0' 
and 'shock' denote quantities before and after the 
shock. The Mach number of the shock, M = v/c s , 



where c s is the velocity of the sound in the pre- 
shock and v is the velocity of the gas with re- 
spect to the plane of the shock (e.g. Landau & 
Lifshitz 1959). Using the above density jumps, 
the Mach number is M = 1.33±g;J| and 1.53±g |J, 
respectively for the east and south-west surface 
brightness discontinuities. From these density 
jumps, we predict temperature jumps T s hock/?o 
= 1.32±g;fl and 1.52±J;||, respectively for the 
east and south-west wedges, which are consistent, 
within error bars, with our measured tempera- 
tures. The pressure is discontinuous at the shock 
front as shown in Figures [5] and [H panel (d) for 
the east and southwest shocks, respectively. 
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4.. 2.1. Merging and 'Sloshing' 

'Sloshing' of the dense cluster gas core is the 
term used to describe the non-hydrostatic gas mo- 
tions induced by a recent merger. When a cluster 
undergoes a perturbation to its gravitational po- 
tential caused by another infalling group or clus- 
ter, the gas core lags behind the cluster the po- 
tential minimum, as they both move toward the 
perturbing object. As the gas core falls back onto 
the potential minimum, it overshoots it and begins 
to oscillate. With each oscillation, the gas core is 
moving against its own trailing gas, producing an 
"edge" in the X-ray brightness which expands out 
from the cluster. This sequence of events is de- 



scribe d in more detail inlAscasibar &, Markevitch 
(|2006t ) and iMarkevitch et al. I <|200lh - The con- 
tinued oscillation of the core gas about the po- 
tential minimum produces a succession of radi- 
ally propagating cold fronts, manifested as con- 
centric edges in the surface brightness distribu- 
tion. These fronts may form a spiral structure 
when the sloshing direction is near the plane of 
the sky and the merger has a non-zero an gular 



momentum ( Ascasibar fc Markevitch 20061 ) 



The observational features of a cold front are 
a large temperature jump across the discontinu- 
ity (the colder gas is closer to the nucleus) with 
no discontinuity in pressure. The gas density, and 
therefore surface brightness, also generally forms 
a spiral pattern. If we assume that the pressure is 
continuous across the discontinuity, the tempera- 
ture jumps should be by factors of 0.67 



013 and 

0-57^013 respectively, based on the derived den- 
sity jumps for the eastern and the south-western 
surface brightness discontinuity. These jumps are 
well within the uncertainties shown in Figures [5] 
and [5J If we assume that we are viewing an 
advanced stage of gas sloshing, the two discon- 
tinuities possibly form a clockwise spiral pattern 
around the core of the cluster because the south- 
western discontinuity is father than the eastern 
discontinuity from the core. In much better ex- 
posed Chandra observations of 'sloshing', however, 
the discontinuity is most prominent in one small 
sector, not equally visible in tw o nearly opposite 



directions ( Johnson et al. 2010) 



To our knowledge, there has been no compre- 
hensive study of individual galaxies of this cluster 
or measurement of their velocity dispersion. The 



nearest galaxy to the cluster center detected in 
the Sloan Survey, SDSS J133850.88+385216.0 (z 
— 0.2442), at a projected distance of —0.3 Mpc 
(1.131') (SDSS DR7:1 lAbazaiian et al. f2009h . has 
no gas associated with it in the Chandra image. 
Additionally, there is no clear pattern of spiral 
structure in the core, only an offset in the positions 
of the discontinuities. However, the relatively low 
quality of the data makes a quantitative statement 
impossible. We conclude that the only compelling 
argument to support the sloshing hypothesis is the 
vague similarity of the discontinuities seen in 3C 
288 with other 'sloshing' systems, and as described 
above there are several lines of evidence to suggest 
they are shocks. A much deeper Chandra obser- 
vation and/or a systematic study of the member 
galaxies of this cluster would provide a definitive 
answer. 

4.3. Nuclear Outburst 

If the surface brightness discontinuities in 
3C 288 are shocks, we can estimate (i) the to- 
tal energy and the age of the outburst, (ii) the 
amount of mass accreted by the central SMBH, 
and (iii) make a strong statement about the ef- 
fects of outbursts on suppression of the formation 
of large amounts of cool gas. The total thermal 
energy of the gas within the core radius of the 
beta-model is —2.0 x 10 60 ergs. This approach 
underestimates the true energy because the shock 
front extends considerably further to the east and 
southwest, implying a faster, stronger shock, en- 
compassing a greater volume in those directions. 
However the value of 2.0 x 10 60 ergs provides a 
conservative estimate of the thermal energy of the 
gas. We estimate the mechanical energy of the 
outburst by two methods. First, assuming that 
the bubbles inflated adiabatically, the total bub- 
ble enthalpy is -9.0 x 10 59 ergs. Only 25% of 
this (— 2.3xl0 59 ergs) has gone to heat the gas. 
We have assumed an enthalpy of 4pV for each lobe 
(i.e. 7=4/3), and that the lobes ca n be modeled as 



cylin ders in the plane of the sky (Raff ertv et al. 
2006). This is the minimum energy of inflation. 
If the inflation were in fact supersonic, the energy 
imparted to the gas could be considerably larger. 
The minimum mean mechanical power of the jet 
is then Pj Ct = 6.1 x 10 44 ergs s , if we assume 
the bubbles ar e buoyantly evolving w ith a mean 
speed of 0.5c s (jChurazov et al. 11200 ll ). 
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We next compute the energy of the outburst 
assuming that the shock is caused by an isotropic 
point explosion. We model the shock as a one- 
dimensional point release of energy into a /3-model 
atmosphere. The parameters of the model, includ- 
ing the energy and age of the burst, are adjusted 
to fit the observed surface brightness profile. This 
model has been used to constrain the shock param- 
eters for NGC 463 6 and other nuclear outbursts 
(|Baldi et al. I [2009). First, for the eastern shock, 



we fit the shock model to the surface brightness 
profile, finding that the shock energy is 4.7 x 10 59 
ergs and the shock age is ~3 x 10 7 yr. The latter 
is better determined, since it depends largely on 
the shock radius and its current speed. Second, 
for the south-west shock, a similar analysis yields 
a shock energy of 7.7 x 10 59 ergs and shock age of 
^4 x 10 7 yr. The two ages are fairly similar and 
this suggests that there is some asymmetry in the 
pressure profile. It is conceivable that the shocks 
arise from two different outbursts. If the shocks 
are part of a front produced by a single outburst, 
the total energy would lie (roughly) between these 
two values, or ^6.2xl0 59 ergs, and similarly the 
shock age would be ^3.5 x 10 7 yr. 

The mass of the central SMBH can be esti- 
mat ed from the K-z da ta for the 3CRR sam- 
ple dWillott et al. I l2003h and is ~4 . x 10 9 



Mr. 



(Hardcast le. Evans, fc Croston 1 2007; M arconi fe Hunt 



2003). Assuming the outburst was powered by the 
gravitational binding energy released by accretion, 
and adopting a mass-energy conversion efficiency 
e = 0.1 and a total outburst energy of between 
-0.9xl0 60 ergs (= ApV) and 1.3 x 10 60 ergs (= 
ZpV + shock-outburst-energy), we find that under 
these assumptions the black hole grew by 



AAf BH = 



(l-e)E 



10 7 Ma 



Here, AMbh accounts for the lost binding energy, 
E is the total energy output in mechanical and 
radiative forms. We ignore the radiation because 
it accounts for a negligible fraction of the current 
power output. This growth in mass corresponds 
to an average growth rate of 0.3 Af© yr -1 over 
the past ~3.5 x 10 7 yr. Thus the current outburst 
is a small (<1%) contribution to the mass of the 
central SMBH. 

Can Bondi accretion of the ICM account for the 



2006)? It is in principle straightforward to regu- 
late in the context of feedback models and cooling 
flows ( Nulse n fc Fabian | l2000t IChurazov et al. 



20021: ISiiacki et al. 1 120071: ISomerville etal. 



200£ 



McNamara et al. I 120091 ) and the X-ray atmo- 
sphere provides a steady supply of fuel. In rel- 
atively low power radio galaxies hosted by gi- 
ant ellipticals, Bondi accretion has been shown 
to be energetically feasible in the sense that 
hot atmospheres probably have a sufficient gas 
density to supply the mass re quired to account 
for the observed j et po wers (IDi Matteo et al. 



20001: lAllen et al. I [20ol iRafferty et al. I 120061 



Hard castle. Evans, fe Croston 1120091 : iMcNamara et al. 
2009h . 

The average gas density and temperature in the 
inner 3.2" (12 kpc) of 3C 288's hot halo is n e = 2.3 
x 10~ 2 cm 3 and 2.9+q'2 keV, respectively. Using 
the black hole mass of Mbh — 4.0 x 10 9 M (see 
above), we find a Bondi accretion rate of 



Mb = 0.012 



0.13 



kT 
2J5 



-3/2 



Mi 



BH 



10 9 



= 2.7 x 1O~ 2 M yr' 1 . 

This value lies roughly an order of magnitude be- 
low the M = 0.3 Mq yr -1 required to power the 
current outburst and taken at face value suggests 
that the current outburst cannot be powered by 
Bondi accretion of hot cluster medium. We cau- 
tion however that this result relies on an extrapo- 
lation of the temperature and density profile into 
the core. Chandra's resolution is not sufficient to 
probe the gas on scales of kiloparsec, and even a 
modest increase in density and decrease in tem- 
perature of the gas on these spatial scales could 
easily balance the Bondi accretion rate with the 
mechanical power of the outburst. For example, 
the gas in the BCGs of hot, non-cool core clusters 
such as Coma i s denser and cooler than that of the 
ambient ICM (jVikhlinin et al. Il2001bl ). The fact 



mechanical power of the outburst ([Allen et al. 



that the shock is detached from the lobes, that 
no X-ray emission is detectable from the central 
AGN, and that the radio galaxy may be transi- 
tioning from FRII to FRI all support the idea 
that inflation of the lobes has slowed perhaps due 
to the energy supply to the jet being greatly re- 
duced or cut off. Therefore, it is not surprising 
that energy released due to the Bondi accretion is 
currently far less than the mechanical power of the 
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outburst, if the AGN power has recently dropped 
significantly. Alternatively, it is possible that the 
current outburst was fueled by accretion of cold 
gas from, for example, a dusty disk. Such accre- 
tion is commonly seen in nearby FRIIs such as 



3C33 (jKraft et al. Il2007bl) . There is no evidence 



of a d usty disk in the HST image (jde Koff et al. 



19961 ) . so it is not clear where this cold gas would 
originate. Hence, in short, it is plausible for the 
Bondi accretion to power the current outburst and 
cold gas could come from a minor merger, but the 
existing data is of insufficient quality to make a 
definitive statement. 

The estimated energy of the shock is roughly 
twice the value of the minimum pV work done by 
the inflation of the lobes, thus demonstrating that 
supersonic inflation of the lobes can play a key role 
in the energy balance of cool core clusters. The 
shock energy is also a significant fraction (~30%) 
of the total thermal energy of the gas within the 
radius of the shock. This possibly suggests that 
the temperature and thermal energy of the gas in 
the core prior to the inflation of the radio lobes 
was at least 30% lower. We conclude that we are 
most likely witnessing AGN feedback in action, 
and that the outburst, which may have been fueled 
by Bondi accretion of cooling gas at the cluster 
center, has likely puffed up the cluster cool core 
to offset radiative losses. 

4.4. Internal Pressure of the Radio Lobes 

Finally, we use the measured pressure profile 
of the gas to determine whether the lobes are at 
or near equipartition. To estimate the equiparti- 
tion magnetic field strength, B eq , in the lobes, we 
make the conventional assumptions that all rele- 
vant features are cylinders with depths equal to 
their radius on the plane of the sky, that the radio 
spectra are power laws from 10 MHz to 100 GHz, 
that the filling factor of the emission is unity, and 
that equal energies reside in the heavy particles 
and the relativistic electrons. With these assump- 
tions, B en is ~50 fj,G o ver the lobes, consistent 



with lBridle et al. I (|1989h . and hence, the equipar- 
tition pressure is ^7.7 x 10~ n dyn cm -2 . Using 
the best-fit models of the /3 profile, we estimate the 
thermal gas pressure at the approximate position 
of the lobes to be ~1.1 x 10~ 10 dyn cm~ 2 . Thus 
for the case of this FRI/FRII transitional object, 
the equipartition pressure of the lobe is roughly 



equal to that of the ambient pressure given the 
uncertainties of roughly 50%. The result does not 
change even if we assume that lobes of 3C 288 do 
not contain an energetically dominant proton pop- 
ulation. This suggests that this 'transitional' ob- 
ject is, at least in this regard, more similar to the 
FRII radio galaxies than the FRIs. It is typically 
found that the equipartition pressures of lobes in 
FRI radio galaxies are orders of magnitude less 
than t he ambient gas, while the P eq 
FRIIs (|Croston et al. Il2005l l2008h . 



P gas for 



5. Conclusions 

3C 288 and its gaseous environment provide a 
laboratory at moderate redshift for investigating 
the interaction between an outburst from a SMBH 
and the surrounding cluster medium. Using the 
42 ks Chandra observations of 3C288, we deduce 
the following: 

1 . We detect two surface brightness discontinu- 
ities in the gas at projected distances of 53.3 
(eastern) and 91.1 kpc (southwestern) from 
the nucleus, which we attribute to shocks 
from the supersonic inflation of radio lobes. 

2. Under the assumption that the disconti- 
nuities are shocks, the gas density jumps 

(/3shock//00 ~ 1.48t°:24 and /°shock//00 « 

1.75^39, respectively with 90% uncer- 
tainties for the eastern and south-western 
shocks), yield shock Mach numbers, 1.33i^i2 
and 1. 531q 21, respectively for the eastern 
and south-western shocks, characteristic of 
a classical shock in a gas with 7 = 5/3. The 
data are not of sufficient quality to detect 
the expected jump in temperature at the 
discontinuity (T shock /T ). 

3. We measure the energy and age of the shocks 
to be ~1.6xl0 60 ergs and 3.5 xlO 7 yrs, re- 
spectively. 

4. The radio lobes are not far from equiparti- 
tion. 

Chandra has detected shocks from the super- 
sonic inflation of radio lobes in nearly two dozen 
galaxies, groups, and clusters, but the discovery of 
shocks reported here is the most distant reported 
to date. It is now clear that feedback between 
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cooling gas and the central SMBH plays a critical 
role in the evolution of early-type galaxies and the 
central regions of groups and clusters. These out- 
bursts probably suppress star formation in massive 
galaxies and are the origin of the exponential de- 
cay in the galaxy mass function at large masses. 
A well-selected Chandra survey could detect a sig- 
nificant number of examples of such phenomena in 
massive systems out to redshift ~0.5. The look- 
back time to 3C 288 is only about 15% of the Hub- 
ble time, so even this observation has not yet be- 
gun to directly study the role of shock-heating in 
the epoch of cluster formation (z ~ 1 and beyond). 
Given the relative faintness and rarity of even the 
most massive clusters beyond z ^0.5, it would be 
difficult to make a detailed study with Chandra 
at and beyond the redshift at which clusters are 
forming. 
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